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Abstract 

The critical properties of the abehan Polyakov loop and the Polyakov loop 
in terms of Dirac string are studied in finite temperature abelian projected 
SU{2) QCD. We evaluate the critical point and the critical exponents from 
each Polyakov loop in the maximally abelian gauge using the finite-size scal- 
ing analysis. Abelian dominance in this case is proved quantitatively. The 
critical point of each abelian Polyakov loop is equal to that of the non-abelian 
Polyakov loop within the statistical errors. Also, the critical exponents are in 
good agreement with those from non-abelian Polyakov loops. 
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I. INTRODUCTION 



Abelian projected QCD has been studied extensively in recent years, for elucidating the 
mechanism of quark confinement |l|J^. The abelian projection of QCD ^ is to perform a 
partial gauge-fixing such that the maximal abelian torus group remains unbroken. Abelian 
monopoles appear as a topological quantity in such a partial gauge fixing, so that QCD can be 
regarded as an abelian theory with electric charges and monopoles. 't Hooft conjectured that 
if the monopoles made Bose condensation, quarks could be confined due to dual Meissner 
effect §. 

There are some evidences on lattices that the abelian theory in the maximally abelian 
(MA) gauge |Q well represents the long range properties of QCD: 

1. Abelian Wilson loops composed of abelian link fields alone can reproduce the full 
{SU{2)) value of the string tension. Furthermore, the abelian Wilson loops written in 
terms of monopole currents also reproduce the value pHlOf • 



2. Polyakov loops written in terms of abelian fields and also in terms of Dirac strings 
of monopoles (monopole Polyakov loops) can reproduce the behavior of non-abelian 
Polyakov loops W^. 



3. A monopole effective action can be calculated. The argument of the energy and the 
entropy indicates that QCD is in the monopole condensed phase 0,|12 . 



These facts are usually called abelian (monopole) dominance in quark confinement and 
suggest that 't Hooft's conjecture is realized in MA gauge. 

Figure |l| shows the non-abelian, the abelian and the monopole Polyakov loops versus (3 on 
24^ X 4 SU{2) lattice The abelian and the monopole Polyakov loops change drastically 



around the critical point (3=2.29 determined from the non-abelian Polyakov loops. The 
abelian and the monopole Polyakov loops appear to be good order parameters. However, 
those curves seem to have different slopes. Their absolute values in the deconfinement phase 
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are also different. Actually, those three, the non-abelian, the abelian and the monopole 
Polyakov loops are quite different operators. 

The critical property of 4- dimensional SU (2) lattice gauge theory is shown to be universal 
to that of 3-dimensional Z2 theory |T^. It is interesting to study what exponents are 
calculated from the abelian and the monopole Polyakov loops at each critical point, since 
Z2 symmetry is not so directly understood in the framwork of monopole dynamics [Q and 
there is no reason the exponents of those different Polyakov loops agree with each other. 
Such a study helps us also to test how good the abelian dominance is quantitatively. It is 
the aim of this work. 



II. DEFINITION OF ABELIAN AND MONOPOLE POLYAKOV LOOPS 



A non-abelian Polyakov loop in SU (2) lattice gauge theory is written in the form 



^ t 



(1) 



where U^{x^t) are SU{2) link variables at space x and at time t. 

After abelian projection is over, we can define abelian Polyakov loops written in terms 
of abelian link variables. The abelian link variables can be separated from gauge-fixed link 
variables 



Uf,{x,t) = C^{x,t)u^{x,t), 



(2) 



where Uf^{x, t) is a gauge-fixed link variable. u^{x, t) is a diagonal matrix composed of phase 
factors of the diagonal components of Ufj_{x,t). We can define an abelian Polyakov loop 



Pabei{xo) = Re[exp(i^6'4(x,t)J4(x,t))]. 



(3) 



x,t 



Here J^{x,t) = 5^,45x,xo and 6'^(x, t) are the angle variables of u^{x,t): 



u^{x,t) 



exp{i9^{x,t)) 







exp(-i6'^(x,t)) ^ 
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(4) 



The abelian Polyakov loop can be decomposed into two parts: a monopole part and a 
photon part |Tl|. An abehan field strength can be written as 



f^,u{x, t) = d^9u{x, t) - du9^{x, t), (5) 
where d^, is a forward derivative. Rewriting this equation, we get 

e^{x, t) = -J2D{x~ x', t - t')dlMx', t') -J2D{x- x', t - t')d^d'Mx', t'), (6) 

x',t' x',t' 

where d'^ is a backward derivative and D{x, t) is a lattice Coulomb propagator which satisfies 
d'^d^D{x,t) = —6xfi6tfi. Then the abelian Polyakov loop (Eq.(^) can be written in terms 
of the abelian field strength: 

Pabeiixo) =Re[eM-i E D{x - x' ,t - t')d'J,^{x' ,t')J4{x,t))]. (7) 

x,t,x',t' 

Here the second term of Eq.(D vanishes owing to d^J^lXjt) = 0. The abelian field strength 
can be separated into two parts: 

fuf, = fuf, + '2iTn^^, (8) 
where riuf^ is an integer and fu^ G [— 7r,7r). Then, rewriting Eq.(^, we get 

Pabei{xo) = Re [exp(-z J2 -^(^ -x',t- t')d'JuA{x', t')Ji{x, t)) (9) 
xexp(-27ri ^ D{x - x\t - t')d'^n^i{x' ,t')Ji{x,t))] (10) 

x,t,x',t' 

= Re[Pi(xo)-P2(xo)]. (11) 

The monopole Polyakov loop, Pmono{.XQ) = RePi(a;o) is composed of Dirac strings of 
monopoles. Pphoton{xo) = ReP2(a^o) only contains the contributions from photons. Suzuki 
et al. UTTH have indicated that 



1. Pabel{Xo) ~ Pmono{Xo) X Pphotonixo) in MA gaUge. 

2. Pabei{xo), Pmonoixo) and Psu{2){xq) vauish for (3 < (3c. 

3. Pphoton is finite from /5 =2.1 to 2.5 and does not change drastically around the critical 
point. 



III. FINITE-SIZE SCALING THEORY 



We calculated the critical exponent of the non-abelian, the abelian and the monopole 
Polyakov loops from a finite-size scaling theory. The singular part of the free energy density 
on X Nt lattice has the following form: 

fs{x,h,N,) = N^'QsixNl'^hNl^^^^'^g.Nf), (12) 

where x = {T — Tc)/Tc. Here the action contains the term hNfL (L denotes the magnetiza- 
tion) and Qi are the irrelevant fields with exponent yi. By differentiating fs with respect to 
/i at /i = 0, we get 

{L){x,N,) = N;^/''QL{xNy'',g,Ny^), (13) 
N,) = N]/''Q^{xNl/\ g,Nf), (14) 
gr{x,N,)=QgXxNl'\9iNf)^ (15) 

where (L), x ^-iid Qr are order parameter, susceptibility and 4-th cumulant, respectively: 

^ = 4En^)^ (16) 

S X 

X = N!m-{L)% (17) 
/ r4\ 

3r = J^- 3. (18) 
Expanding these equations with respect to x, we have at x = 

{L){x = 0, Ns) = iV7'3/^(co + caiV,-), (19) 
X{x = 0, iV,) = iV;/-(co + csiVr ), (20) 
gr{x = 0,Ns) = co + csN;'^, (21) 

where we take only the largest irrelevant exponents {—uj) into account. We can calculate the 
critical point from the fit to the A'^g-behavior of those equations. The critical point can be 
defined as the point where a fit to the leading A^'^-behavior has the least 10 • Actually, 
the leading A'^g-behavior of Eq.(|l9D and of Eq. (|20|) is given by 



\n{L){x = 0,Ns) = -^\nNs + \nco, (22) 
lnx(x = 0,Ar,) = ^InA^, + lnco. (23) 

From the fits to Eqs.(|2l|), (|2|) and (H), we can find the position of the critical point /?c, 



and obtain the values of /5/z/, 7/z/ and at Pc simultaneously. Here g"^ is the value of gr 
on the infinite volume and is denoted by Cq in Eg . (pT]) . We also considered the derivatives 
of the observables with respect to x. The A^^-behavior of each derivative is given by 

^{x, N,) = Ni^-^y^Q'^{xNl/^), (24) 



N,) = Ni^+'<^'^Q'^{xNl/^), (25) 



^(x,iV.)=iVy^Q;,(^Ary^)- (26) 
The leading A^^-behavior of each equation at the critical point is 

In ^ (x = 0, iV,) = 1^ In + In ^o, (27) 
ox V 

In |^(x = 0, iV,) = IniV, + In do, (28) 

ax V 

ln^(x = 0,A/'J = -InAT, + lnrfo- (29) 
ax V 

Hence, (1 — /?)/z/,(l + 7)/// and 1/z/ can also be evaluated from the fits to those equations. 



IV. RESULTS AND DISCUSSIONS 

We performed numerical calculations on A^^ x 4 lattices, where A^^ =8, 12, 16 and 24. 
The standard SU{2) Wilson action was adopted and abelian link valuables were defined in 
MA gauge. We calculated the observables 

L = ^sT.Pix), (30) 

where P{x) denotes Psu(2){x), Pabei{x) and Pmono{x). Actually, Xv = Ng{L'^) was calculated 
instead of X; because Xv is equal to x below the critical point |T^. The values of the 
observables at various (3 are needed in order to calculate the derivatives with respect to 
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X, where x = {P — Pc)/Pc- We used the density of state method(DSM) |[T6|JT7[] . First we 
performed Monte-Carlo simulations at P = Pq, and then calculated the following averages: 

N{s) 

0{s) = E Ok/N{s), (31) 
fc=i 

where s is the value of the action, N{s) the number of the configurations whose action 

has the same value of s, and Ok the observable obtained from the k-th configuration. The 

expectation value of the observables in the vicinity of Pq is then given by 

^ j:sN{s)0{s)exp{-Ps + Pos) 

EsN{s)exp{-Ps + Pos) ' ^^^^ 

where Pq =2.2988 was adopted and P €[2.2980,2.3000]. (L), x and Qr at Pq were calculated 

every 50 sweeps after 2000 thermalization sweeps. The number of samples was 100000, 

except on 24^ x 4 lattice (47000 in the case). The errors were determined according to the 

Jackknife method dividing the entire sample into 10 blocks (4 blocks on 243 X 4 lattice). 

We estimated the critical point Pc from the method. The data of our DSM results 



were fitted to Eqs.(pi])-(^) and Eqs.(p7D-(p9|) at each p. The number of input data was 



2 and that of fit parameters was 2 {00 in Eq.(21) was fixed to 1 in accordance with Engels 



et al. [|3). Figure || describes the typical curves of /Nf versus p. Here the number of 
degrees of freedom, Nj is 2. Each curve in Fig. |^ is smooth and has its minimum value. 
Table | shows the positions of minimal l^f fo^' observables obtained. Almost all the 
X^in/Nf are small. However, the xlninl-^f were not seen from our dgr/dx fits. Furthermore, 
our two-parameter fits were not so good in the cases of (L), x, d{L)/dx and dx/dx from the 
monopole Polyakov loops. Then we used Eq.(|r^) and Eq.(pil|) for their fits which contained 
three parameters to be fitted. The values of a; in Eq. (|TP|) and Eq. ([2^) were chosen in such 
a way that the values of Xmin/^f became as small as possible. 
Averaging the obtained minimal positions of Xmin/-^f^ 

psu{2) ^ 2.29940(20), (33) 
^abei ^ 2.99962(26), (34) 
/?™ = 2.29971(23). (35) 



The critical points obtained from the abehan and the monopole Polyakov loop are very close 
to the non-abelian critical point as expected from Fig. 

Table [I^ lists the critical exponents on each critical point in the non-abelian, the abelian 
and the monopole case. The non-abelian exponents obtained by Engels et al. [|T5|, the 
exponents of 3-dimensional Ising model and those of f/(l) |T^ are also shown. The errors 



were caused by fluctuations of the interpolated DSM data and by uncertainty of each Pc- See 
also Fig. Those critical exponents seem to be reliable because of the following reasons: 
three u's obtained from three different fits are within the statistical errors; hyperscaling 
relations are well satisfied; non-abelian exponents obtained are consistent with those of 
Engels et al. [1^ . Table || shows the following notable results: 



1. The critical exponents in the abelian and the monopole case are in agreement with 
non-abelian exponents within the statistical error. 

2. Those critical exponents agree with those of Z2 rather than those of U{1). 

The abelian (monopole) dominance in quark confinement is proved quantitatively in this 
case. 

There remain some problems to be studied further. The critical points obtained are 
outside the one-sigma error bar of Engels et al. |jl5|; minimal /Nf were not seen for some 



fits; some values of Xmin/^f ^(1) which are not small enough. These problems seem to 
reflect a lack of statistics. More samples may be needed especially on 24^ x 4 lattice. 

The simulations of this work were carried out on VPP500 at Institute of Physical 
and Chemical Research (RIKEN) and at National Laboratory for High Energy Physics 
at Tsukuba (KEK). This work is financially supported by JSPS Grant-in Aid for Scientific 
Research (B) (No.06452028). 
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TABLES 

















\L„M7 


9r 


SU(2) 


2.29952 


1.49 


dgr/dx 


SU(2) 


- 


< 0(10-^) 




abel 


2.29936 


1.27 




abel 


- 


< o(io-i) 




mono 


2.29971 


0.17 




mono 


- 


< 0(10-^) 


(L) 


SU(2) 


2.29960 


1.33 


d{L)/dx 


SU(2) 


2.29920 


0.004 




abel 


2.29984 


0.95 




abel 


2.29938 


0.017 




mono 




< 0(1) 




mono 


2.29918 


G X IQ-'^ 


X 


SU(2) 


2.29946 


1.33 


dx/dx 


SU(2) 


2.29924 


0.005 




abel 


2.29986 


0.75 




abel 


2.29964 


0.061 




mono 




<0(1) 




mono 


2.29992 


6 X 10-6 



TABLE I. The positions of minimal x^/^f the value of Xmin/-^f non-abelian,the 
abelian and the monopole cases. 
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SUi2) 


abel 


mono 


Engels et al. |15 


Ising PI 

II 11 


U(l) pi 


PI 




n 485('99'l 


598('64') 


595('8'l 


n ^i8i'7'i 




(1 


1.117(27) 


1.138(10) 


1.091(84) 


1.085(14) 


1.072(7) 




V 


0.617(16) 


0.616(12) 


0.617(57) 


0.621(6) 


0.6289(8) 


0.67 


(i 


0.311(19) 


0.299(19) 


0.326(69) 


0.326(8) 


0.3258(44) 


0.35 




1 Q77('9Q') 


9 n95('34'l 

^ .\JZjO\*J'-t J 


1 QQ1 (RR) 


1 Qdd(^ 3") 

± . c/'-t'-ty -Lkj 1 


1 970^ 1 




\^ ^ 1)1 


3 6n0f38') 


R 646 ('44) 


3 608 ('93') 


3 555fl5') 


3 56n('ll'l 




V 


0.616(25) 


0.617(29) 


0.618(68) 


0.621(8) 


0.6289(8) 




7 


1.218(68) 


1.249(81) 


1.23(19) 


1.207(24) 


1.239(7) 


1.32 


7/1/ + 2/3/1/ 


2.985(47) 


2.995(56) 


3.05(15) 


2.994(21) 


3.006(18) 






1.447(41) 


1.438(42) 


1.438(41) 


1.403(16) 


1.41 






0.633(13) 


0.621(14) 


0.600(13) 


0.630(11) 


0.6289(8) 





TABLE II. The critical exponents calculated from the non-abelian, the abelian and the 
monopole Polyakov loops at each critical point. 
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FIG. 1. Non-abelian, abelian and monopole Polyakov loops in SU{2) QCD on 24^ x 4 lattices. 
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FIG. 2. X^/^f from g,. fits and from d{L)/dx fits versus /3 in the non-abelian, the abehan and 
the monopole cases. The number of degrees of freedom, Nf is 2. 
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Critical exponents 
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FIG. 3. Critical exponents of non-abelian, abelian and monopole Polyakov loops in SU{2) 
QCD. 
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